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In the search for the nature of dark energy most cosmological probes measure simple functions 
of the expansion rate. While powerful, these all involve roughly the same dependence on the dark 
energy equation of state parameters, with anticorrelation between its present value wo and time 
variation Wa- Quantities that have instead positive correlation and so a sensitivity direction largely 
orthogonal to, e.g., distance probes offer the hope of achieving tight constraints through complemen- 
tarity. Such quantities are found in strong gravitational lensing observations of image separations 
and time delays. While degeneracy between cosmological parameters prevents full complementarity, 
strong lensing measurements to 1% accuracy can improve equation of state characterization by 15- 
50%. Next generation surveys should provide data on roughly 10^ lens systems, though systematic 
errors will remain challenging. 



I. INTRODUCTION 

Dark energy poses a fundamental challenge to our un- 
derstanding of the universe. The acceleration of the 
cosmic expansion discovered through the Type la su- 
pernovae distance-redshift relation can be interpreted in 
terms of a new component of the energy density possess- 
ing a substantially negative pressure. Further observa- 
tions from the cosmic microwave background radiation 
indicate the spatial geometry of the universe is flat; in 
combination with the supernova measurements this im- 
plies that the unknown "dark energy" comprises roughly 
70% of the total density, in concordance with large scale 
structure data indicating that matter contributes approx- 
imately 30% of the critical density. 

Such a weight of dark energy leads to its dominance 
of the expansion, causing acceleration, restricting the 
growth of large scale structure, and holding the key to the 
fate of the universe. But apart from its rough magnitude, 
its nature is almost unknown - whether it arises from the 
physics of the high energy vacuum, a scalar field, extra 
dimensional or "beyond Einstein" gravitational effects, 
etc. One way to gain clues to the underlying physical 
mechanism is to characterize the behavior of dark energy 
in terms of its equation of state ratio as a function of red- 
shift (essentially a time parameter), This is con- 
ventionally interpreted in terms of the ratio of the dark 
energy pressure to energy density, but can also be used 
as an effective parameterization to treat generalizations 
of the cosmology framework of the Friedmann equations 
of expansion 1] or simply in terms of the expansion rate 
and acceleration itself j^]. 

To obtain the clearest focus on the class of physics re- 
sponsible for the accelerating universe, we seek the tight- 
est constraints on the equation of state. This must come 
not merely from high statistical precision of a probe, 
but from robust control of systematic uncertainties. The 
greatest accuracy and confidence in the measurements 
will come from independent crosschecks and complemen- 
tarity between methods of probing the cosmology. Many 
studies have considered such complementarity between 
probes (e.g. 0,1111111111), with a promising future for 



next generation surveys carrying out such measurements. 
However, all the methods of supernova distances, cos- 
mic microwave background power spectrum, weak grav- 
itational lensing, cluster counts, and baryon oscillations 
possess a similar fundamental dependence on the equa- 
tion of state through the Hubble parameter, or expansion 
rate. 

This paper investigates whether a truly complemen- 
tary probe exists that has nearly orthogonal dependence 
to the previous ones in the plane of the equation of state 
parameters of value today, wq, vs. time variation, Wa- 
Such a probe, if practical, would offer a valuable contribu- 
tion to uncovering fundamental physics and deserve fur- 
ther consideration among next generation experiments. 
Section ^1 considers the characteristics such a method 
would possess and identifies two promising candidates re- 
lated to strong gravitational lensing. In ^IHII we analyze 
the sensitivity of these probes to the cosmological pa- 
rameters and the constraints and complementarity they 
offer. In the conclusion we summarize the prospects for 
strong lensing as a cosmological probe and discuss some 
issues regarding surveys and systematic uncertainties. 



II. COMPLEMENTARITY IN THE EQUATION 
OF STATE PLANE 

Astronomical observations involving distances and vol- 
umes all follow from the metric in a simple, kinematic 
way ^9]. For convenience, the dependence can be written 
in terms of the conformal distance interval 

drj = dt/a{t) = -dz/H{z), (1) 

where a{t) = (1 -I- z)~^ is the scale factor as a function of 
proper time, equivalently parametrized in terms of red- 
shift z, and H = a/a is the Hubble parameter. So in 
observing some source at redshift z the cosmological in- 
formation is carried by the quantity H . 

Varieties of distance observations - angular diameter, 
luminosity, proper motion - all contain the same infor- 
mation. One exception is the parallax distance (cf. [lo| . 
§3.2) which involves not only H but the spatial curva- 
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ture k as well; however this is not practical for cosmol- 
ogy and CMB measurements strongly indicate a flat uni- 
verse. Another subtlety involves exotic models that break 
the thermodynamic, or reciprocity, relation between the 
various distances [TH IT^ by violation of Liouvillc's 
theorem, but this arises from photon properties and not 
cosmology. 

Volume elements, and hence numbers of sources, are 
built up out of distances and so similarly involve the 
Hubble parameter. From the Friedmann equation the 
relation between the expansion rate and the dark energy 
equation of state is 



flDEe 
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where Hq is the present value of the Hubble parameter 
(the Hubble constant), and for a flat universe Ude — 
1 — ^Im- As discussed in this can be written more 
generally as defining an effective equation of state 
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where 5H^ encodes our ignorance of the right hand side 
of eq. 121 after the first, matter density term. 

Making w{z) more positive (holding ifp and Q.m fixed) 
increases the expansion rate in the past, and so decreases 
the acceleration (since the rate today is fixed). Distances 
will be less, sources appear brighter, etc. If we try to 
characterize the nature of the dark energy in the simplest, 
nontrivial way, through the value of its equation of state 
today and its increasing or decreasing positivity into the 
past, we see that cosmological measurements, depending 
only on iJ, don't really care in a gross sense how w{z) was 
more positive ~ either increasing the present value wq or 
increasing the time variation Wa d'w/d\'a{\ + z) does 
the trick. Thus, for all such probes the quantities wq, Wa 
will be anticorrelated; an increase in one can be (at least 
partially) offset by a decrease in the other. If we plot 
constraints from astronomical data in the wq — Wa plane 
(marginalizing over any other parameters), the probabil- 
ity contours will show a degeneracy direction tilted coun- 
terclockwise from the vertical. This is well known and 
illustrated for a number of different probes in 0|; note 
also that a "bare" H of course has the same dependence 

However, this has the fundamental implication that 
complementarity between methods can only be partial in 
this physically key parameter plane; different data sets 
cannot be orthogonal. This restricts the ability to con- 
strain the equation of state and understand the nature of 
dark energy. What would be valuable is construction of 
a cosmological probe whose sensitivity lies clockwise of 
vertical, into the unexplored half-plane of positive corre- 
lation between wq and Wa- Note that this argument does 
not depend on the specifics of how wq, Wa are defined, 
but for concrctcncss we adopt the parametrization 

w{z) = Wq + Wail ~ a) ^ Wq + WaZ/ (1 -|- z), (4) 



successful in fitting smoothly varying dark energy models 
[l6| . A characteristic time variation of the equation 
of state is given by w' = dw/dln{l + z)\z=i = Wa/^. 

Figure ^ illustrates this anticorrelation for the 
distance- redshift probe (dashed ellipse). A hypotheti- 
cal probe that involves positive correlation between the 
equation of state parameters (dotted ellipse) can give 
nearly orthogonal constraints and so the joint probabil- 
ity contour (solid ellipse) between the two types of probe 
tightly characterizes the equation of state. Furthermore, 
since reduced data quality, e.g. from systematic uncer- 
tainties, often lengthen the confidence contour along the 
degeneracy direction (long dashed ellipse), an orthogonal 
probe immunizes against systematic errors. However, de- 
generacies with other cosmological parameters interfere 
with both these desiderata, as discussed in ^IIII 




FIG. 1: This cartoon illustrates the value of orthogonality 
in the equation of state plane for tightly constraining dark 
energy parameters. The short dashed ellipse represents a 
distance-redshift probe, or any other that has an anticorre- 
lation between the equation of state parameters. Adding a 
positive correlation probe (dotted) gives a joint contour (in- 
ner solid curve) with considerably smaller uncertainties. This 
also helps to reduce the impact of systematics (long dashed 
curve and outer solid curve). However the situation is more 
complicated when marginalizing over a larger set of parame- 
ters. 



The growth of large scale structure, or equivalently the 
evolution of gravitational potentials, and hence, say, the 
number of galaxy clusters of a certain mass, might seem 
to offer a different dependence than distances. But the 
growth equation also has H as its predominant ingredient 
and indeed in [3| we see that the growth factor still 
possesses the anticorrelation. An exception might occur 
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for models where the dark energy (or alternative gravity 
theory) is inhomogeneous itself and can act as a source to 
the matter density perturbations, but this is not expected 
to occur for scalar field dark energy on subhorizon scales 

If simple measures of distances or products of distances 
all give anticorrelation, yet they involve differing levels of 
sensitivity. This provides the hope that opposing these 
quantities through ratios might break the pattern of de- 
generacy. Indeed this was found in 0] on consideration 
of the cosmic shear test^ originated by 0|. Here the 
key quantity was H{z)dA{z) ~ H{z) J dz/H{z), where 
d-A is the angular diameter distance. For a certain red- 
shift range, z ~ 1.3 — 2.3 (for = 0.3), the competition 
between the two ingredients causes a positive correlation 
between wo and Wa- A critical pitfall however was that 
systematic uncertainties could blow up the constraint 
contour in a preferred direction such that the resulting 
contour actually exhibited an anticorrelation (compare 
Figs. 5 and 8 of ;4i]). Another indication of positive cor- 
relation was found between discrete values of w{z) in 
neighboring redshift bins by for ratios of distances 
entering weak gravitational lensing calculations. 

Here we investigate distance ratios that appear in 
strong gravitational lensing measurements. Specifically, 
we consider ris/vg, related to the angular separation A0 
between multiple images of a source, ris/{rirs), related 
to the time delay At between the images, and riri^/rs, 
connected with the length scale of caustic properties, or 
the cross section, of the lensing. Here r; is the comoving 
distance to the lens, to the source, and ris between the 
source and lens; in a flat universe ris — rs — ri. 

Figure [21 illustrates the sensitivity of each of these ra- 
tios to the cosmological parameters as a function of red- 
shift. We have taken the redshift to represent the lens 
redshift zi and for convenience fixed the source redshift 
to be Zs — 2zi (since that gives roughly the greatest prob- 
ability of effective lensing). Generally the greatest sensi- 
tivity is to the magnitude of the matter (or dark energy) 
density, then to the present value of the dark energy equa- 
tion of state ratio, and least to the time variation of the 
equation of state. 

Of most interest to our investigation are the crossings 
from negative to positive sensitivity, i.e. at some redshift 
the distance ratio switches from diminishing to growing 
as we change the value of a parameter — {flm, wq, Wa} 
(increasing it, say). Since these crossings occur at differ- 
ent redshifts if the parameter is wq than if it is Wa , then 



A note on the name: this probe involves the observed shearing of 
a sphere due to the properties of spacetime itself. Despite a per- 
fectly isotropic space, a sphere will appear sheared because the 
spacetime is not isotropic - directions along and perpendicular 
to the null geodesies differ. This shear is due to cosmic prop- 
erties, unlike shears in gravitational lensing due to anisotropics 
not in spacetime but in space, though the latter is unfortunately 
sometimes called cosmic shear. 
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FIG. 2: The sensitivity of the distance ratios to the cosmolog- 
ical parameters d — {0,^,1^0,111 a} are encoded in the deriva- 
tives plotted here, normalized to 1% fractional measurement 
of the distance ratio. The larger the absolute magnitude of 
the derivative at a particular redshift, the more constraining 
the observations there. 



the correlation between these quantities can shift from 
negative to positive. We see that positive correlation oc- 
curs for tis/ts for zi « 0.9 — 1.7 and for ris/(rirs) for 
zi ~ — 0.6; zero crossing in equation of state variables 
and hence positive correlation does not occur for riria/ra. 
(The exact values will depend on the fiducial cosmology, 
here taken to be Q.ra = 0.3, wq — —1, Wa — 0.) 

For the ratio rig/r^ we expect rapid evolution in the 
Wo ~ Wa degeneracy direction as the contour will rotate 
from horizontal at z = 0.9 (no sensitivity to wq, so lying 
parallel to the WQ-axis) to vertical at z — 1.7 (no sensi- 
tivity to Wa)- Fig. 131 shows this behavior in a "flower" 
plot, with idealized precision to better illustrate the de- 
generacy direction. A strong lensing survey measuring 
image separations from lenses within this redshift range 
should possess good internal complementarity and the 
ability to constrain the dark energy equation of state. 
Furthermore, in this range strong lensing will have near 
orthogonality with supernova distance and weak gravita- 
tional lensing surveys and so be a valuable complemen- 
tary probe. 

The ratio ris/{rirs) appears somewhat less promising. 
Its sensitivity to Wa is much weaker in the key redshift 
range and the evolution of the correlation between Wq and 
Wa much less dramatic. In particular the contour never 
tilts far from the vertical. This means it is not as com- 
plementary either internally or with other cosmological 
probes. But it does possess three interesting character- 
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FIG. 3: Confidence contours for an idealized experiment mea- 
suring the ratio ris/rs rotate depending on the redshift. De- 
generacy directions follow the Fisher sensitivities from Fig. 
121 with zero crossings of a parameter giving contours paral- 
lel to that axis. Note the strong positive correlation between 
equation of state variables for zi — 1.3, and hence strong 
complementarity with, e.g., distance measurements. 



to the parameter wq a.t z = 0.3 is -6.5 so the best possi- 
ble estimate of Wq from such a single, 1% measurement is 
1/6.5=0.15. That is, the unmarginalized uncertainty is 
cr{wa) =0.15. This will be degraded by degeneracies with 
other parameters, or by worse precision in the measure- 
ment, and improved by additional measurements, subject 
to some systematics floor. 

In fact, degeneracies with additional parameters wash 
out the simple orthogonality of Fig. ^ The parame- 
ter phase space is actually 3-dimensional (from 17™) or 
higher and the joint probabilities do not merely trace 
the intersection of the contours in the equation of state 
plane. This weakening of apparent complementarity also 
means that the magnitude of the sensitivity with redshift 
matters as well as the degree of correlation, and this can 
alter the favored redshift range. All these effects need to 
be taken into account. 

Here we seek to establish what improvements the com- 
plementarity of strong lensing data with other probes 
such as distance-redshift or weak lensing data can realis- 
tically provide. Next generation precision in all three 
of these come in one package: the proposed Super- 
nova/Acceleration Probe (SNAP: With some lO"^ 
galaxies observed in its deep survey (down to AB mag 
R = 30) and 10* in its wide survey (to AB mag R = 27.5), 
the canonical estimate of one strong lens system per thou- 
sand galaxies provides an anticipated wealth of data. 
Other valuable data sets will come from LSST and 
LOFAR m. 



istics: its overall sensitivity to the equation of state is 
reasonably large at z > 1, the positive correlation region 
is at low redshifts where observations are easier (though 
there is less volume, hence fewer lens systems), and it has 
an odd null sensitivity with respect to Qm at z = 0.35. 
This last is where its sensitivity to Wa is maximal in the 
positive correlation region and eases degeneracy with the 
matter density and hence reduces the need for a tight 
external prior on it. So time delay measurements should 
not be wholly neglected as a possibly useful probe. 



III. CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

So far we have only concentrated on the degeneracy di- 
rection. What about the actual magnitude of constraints 
capable of being imposed on the cosmological parame- 
ters? If we consider one parameter at a time, fixing 
the others, we can obtain a lower limit on the parameter 
estimation uncertainty. This is given by 



se > 



I 9 [Distance Ratio] 



de 



(5[Distance Ratio]. (5) 



In Fig. 13 we took a fractional measurement error of 1%. 
Thus, for example, the sensitivity for Hs/r^ with respect 



A. Image Separation 

If we consider an optimistic strong lensing measure- 
ment precision of 1% in a redshift bin of 0.1, whether 
through statistical means, e.g. 100 lens systems of 10% 
precision, or as a systematic floor, we find that the pa- 
rameter uncertainties are large if we do not a priori fix 
some of the variables. Even employing the orthogonal 
combination of measurements of the strong lensing sepa- 
ration variable S = ris/rg at z = 0.9 and z = 1.7 (cf. Fig. 
I^J, the estimations are a{wo) — 0.15, a{wa) = 0.53 even 
with fixed, and 0.28, 0.93 with a prior on of 0.03. 
The lack of sensitivity of the strong lensing separation 
to the dark energy equation of state and, as mentioned 
above, the degeneracy with the other parameter, flm, 
prevent as well its complementarity with other probes 
from being as effective as might have been hoped from 

m 

Figure 0] shows contours in the wq — Wa plane, 
now marginalized over flm, for various combinations of 
probes. This considers 1% measurements at z — 0.9, 1.3, 
1.7. The improvements in parameter estimation remain 
modest: about 15% in wq and Wa when added to SNAP 
supernova distance data, 19% and 3% respectively when 
added to the supernova plus weak lensing data set, 25% 
and 13% when added to supernova plus Planck CMB 
data, and 17% and 6% when added to all three data 
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sets. If strong lensing S measurements were available 
from z — 0.1 — 1.7, say (recall this is the lens redshift, 
with the source assumed to lie at twice this redshift), the 
improvement to the triple data set increases to 30% and 
15%. If the measurement floor for the strong lensing was 
instead at 2%, then the constraints in conjunction with 
the other three probes only gain due to strong lensing 
by 6% and 2% for z = 0.9, 1.3, 1.7 or 17% and 7% for 
z = 0.1 - 1.7. 



provide improvement by 21% and 9% in wq and Wa- 
Adding them to all three other data sets helps by 13% 
and 9%, or 22% and 15% if the strong lensing measure- 
ments extend from z — 0.1 — 1.7. Weakening the precision 
to 2% adjusts this last case to 7% and 5%. Note that no 
external prior on Hq is required; it is determined from 
the data to better than 1%. 



C. Image Separation and Time Delay 
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FIG. 4; Estimates of the dark energy equation of state pa- 
rameters from various combinations of cosmological probes 
(SN=supernovae, WL=weak lensing, CMB=cosmic mi- 
crowave background) improve modestly with the addition of 
strong lensing image separation measurements. Outer con- 
tours of each pair of line types include the labeled data sets 
as will be provided by SNAP and Planck; inner contours add 
the 1% strong lensing constraint. 



Simultaneous addition of both separation and time de- 
lay information allows further improvements. Relative to 
the supernova only data, these are 36% and 41% for 1% 
strong lensing data at z = 0.1, 0.3, . . . 1.7; added to all 
three other probes they yield tighter constraints by 31% 
and 17%. With only 2% precision the improvements in 
the latter case are 15% and 7%. We show the parameter 
constraints in Fig. [S] for the three probes without strong 
lensing, with strong lensing for the fiducial redshift bins, 
for a reduced strong lensing redshift range, and for the 
weaker 2% strong lensing precision. 

The highest redshift bins are not that powerful but 
it is important to include the z = 0.1 bin to constrain 
Hq, which is strongly degenerate with wq and Wa- Even 
with Hq known to better than 1% from the data, the 
degeneracy still plays an important role since the strong 
lensing time delay probe is so much more sensitive to Hq 
than to the equation of state. 

As usual, taking the fiducial cosmology to be that of 
a cosmological constant tends to underestimate both the 
sensitivity and complementarity of the probes For 
comparison Fig. |S1 includes contours for a supergravity 
inspired dark energy model with time varying equation 
of state: S(UGRA) with wo = -0.82, Wa = 0.58. Here 
strong lensing offers more dramatic improvements, by 
49% in estimating wq and 54% in Wa for 1% precision 
and 29% and 22% for 2% precision. Even at 2% preci- 
sion the final uncertainty from the four probes in comple- 
mentarity achieves an impressive constraint on the time 
variation of a{w') = 0.05. 



B. Time Delay 

The other strong lensing quantity we consider is T = 
Hs/ {riVs), related to the time delay between images. This 
introduces an extra parameter - the Hubble constant 
Hq, which we marginalize over. Again, measurements 
of strong lensing by itself provide only poor estimates 
of the cosmological parameters, even for combinations 
of redshifts where the degeneracy directions are comple- 
mentary. The insensitivity to at low redshift noted in 
iO does not help since the overall weak dependence and 
degeneracy of the equation of state variables prevents 
precision constraints. 

In complementarity with supernova distances, strong 
lensing T measurements of 1% precision at z = 0.1 — 0.6 



IV. CONCLUSION 

While strong lensing does not completely fulfill the 
promise of orthogonal constraints on the dark energy 
equation of state that would provide vastly improved 
parameter estimates, it does offer some complementar- 
ity, especially for models with time varying equation of 
state. Furthermore, the data resources required will be 
innate within the deep and wide optical and near infrared 
surveys of SNAP; abundant strong lensing data should 
also come from LOFAR at radio wavelengths and both 
strong and weak lensing from LSST in the optical. We 
find that strong lensing adds reasonable value to dark 
energy constraints when the measurements reach the 1% 
precision level. This is in accord with other calculations 
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FIG. 5: Estimates of the dark energy equation of state param- 
eters from supernovae, weak lensing, and CMB measurements 
can be further tightened by adding observations of both strong 
lensing quantities (SL^): image separations and time delays. 
The outer, long dashed ellipse with no strong lensing and "A" 
in its center corresponds to the outer solid ellipse of Fig. 2] 
Strong lensing data at 1% precision adds valuable comple- 
mentarity (solid ellipse), even over a reduced redshift range 
(dotted). Higher systematic uncertainties on strong lensing 
(2%: short dashed], however, strongly diminish its useful- 
ness. As usual (cf. |g), dark energy models with time varying 
equation of state such as the SUGRA model (ellipses centered 
around "S") show enhanced complementarity. 



(e.g. |2J,|25|)- Precision arises from statistics subject to 
a systematics floor. 

The challenge in the next generation will not be the 
search for sufficient strong lensing data, it will be achiev- 
ing systematics limits sufficiently low to allow the statis- 
tical wealth and the value of strong lensing complemen- 
tarity to come into play. Currently 1% sounds rather 
optimistic for a bound to systematic uncertainties. The 
main contribution is likely to come from our ignorance of 
the lensing mass model. For example, within a singular 
isothermal sphere model 1% measurement of the distance 
ratio S requires knowledge of the velocity dispersion to 
0.5% 26] . Cross correlation of different sources with the 
same lens to remove systematics, a la for weak lens- 
ing, would not work, both because of the rarity of strong 



lensing and because sources at different distances and 
positions would probe different parts of the lensing mass 
distribution. 

One might hope to constrain the lensing mass distri- 
bution by making use of the simultaneous weak lensing 
information that SNAP or LSST would provide, but this 
will be of little use according to jl^ . Basically the weak 
lensing measurements of shears average over a larger 
range of scales (or multipoles) than the convergence con- 
tribution from the mass to the strong lensing. It is like 
trying to locate a pin with thick gloves. Nongaussian 
effects also need to be incorporated in a more rigorous 
treatment of strong lensing. Relaxation of the simplifi- 
cation that sources lie at redshifts Zg = 2zi is unlikely to 
change significantly the results shown here. 

The advantage of image separations and time delays 
is that they are well defined and used as markers of the 
cosmic geometry. One could also use the statistics of 
lensing systems, e.g. the number of lenses out to some 
redshift (e.g. [l^, 113,1^), but then one must deal with 
astrophysical scatter and observational bias and incom- 
pleteness. Furthermore, lensing counts depend on mass 
density growth and volume factors that take us further 
from our original goal of finding orthogonal probes in the 
equation of state plane. 

Since systematics are the limiting factor, we can at- 
tempt to seek out special lensing systems that ameliorate 
these uncertainties, such as strongly lensed calibrated 
candles (e.g. Type la supernovae |33,|33), time evolving, 
differentially amplified sources (e.g. Type II supernovae 
34 13^]), or relatively clean Einstein cross or ring images. 
With the wealth of future data such subsampling may be 
practical. 

While interrelations between cosmological parameters 
entering into distances and the expansion rate pre- 
vent clean orthogonality between cosmological probes, 
the various methods retain sufficient complementarity 
to both crosscheck each other and further tighten con- 
straints on the nature of dark energy. If the systematics 
challenge can be met, in special systems at least, strong 
lensing merits further consideration in joining the cosmo- 
logical toolbox of next generation probes. 
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